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Eclipsing Light Curve for Accretion Flow around Rotating 
Black Hole and Atmospheric Effects of Companion Star 
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ABSTRACT 

We calculate eclipsing light curves for accretion flows around a rotating black hole 
taking into account the atmospheric effects of the companion star. In the cases of no 
atmospheric effects, the light curves contain the information of the black hole spin 
because most of X-ray photons around 1 keV usually comes from the blue-shifted part 
of the accretion flow near the black-hole shadow and the size and the position of the 
black-hole shadow depend on the spin. In these cases, when most of the emission comes 
from the vicinity of the event horizon, the light curves becomes asymmetric at ingress 
and egress. We next investigate the atmospheric absorption and scattering effects of 
the companion stars. By using the solar-type atmospheric model, we have taken into 
account the atmospheric effects of the companion star, such as the photoionization 
by HI and Hel. We found that the eclipsing light curves observed at 1 keV possibly 
contain the information of the black-hole spin. However, in our atmospheric model, 
the effects of the atmosphere are much larger than the effects of the black-hole spin. 
Therefore, even in the case that the light-curves contain the information of the black 
hole spin, it may be difficult to extract the information of the black hole spin if we 
do not have the realistic atmospheric profiles, such as, the temperature, the number 
densities for several elements. Even in such cases, the light curve asymmetries due to 
the rotation of the accretion disc exist. Only when we have the reliable atmospheric 
model, in principle, the information of the strong-gravity regions, such as, the black 
hole spin, can be obtained from the eclipsing light curves. 

Key words: accretion: accretion discs — black hole physics — X-rays: stars. 



1 INTRODUCTION 



Elucidating the nature of the physics of extremely strong 
gravity around black holes is one of the greatest challenges 
in astrophysics in this century. In order to obtain physical in- 
formation of the strong gravity region in the vicinity of black 
holes, direct imaging of a black hole is an extremely promis- 
ing method. The emission from the vicinity of the event hori- 
zon contains information about the physical parameters of a 
black hole. Since the apparent shapes of the black hole shad- 
ows and the appearances of the circular orbit around the 
black hole are deformed by the effects of the black hole's ro- 
tation and charge (e.g. Cunningham & Bardeen 1972, 1973, 
Bardeen 1973, Takahashi 2004, 2005, Zakharov et al. 2005), 
the size and the position of the most luminous parts gen- 
erally existing in the vicinity of the shadow depend on the 
physical parameters of black holes. 

In general, methods for probing strong-field gravity near 
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supermassive black holes or stellar-mass black holes are quite 
limited. Although the past VLBI observations give the im- 
ages of the galactic centers with the spatial resolution within 
100 Schwarzschild radii (e.g., for Sgr A*, Shen et al. 2005, 
for M87, Junor, Biretta & Livio 1999), the direct map- 
ping of the shadows of both supermassive and stellar-mass 
black holes has not been performed so far. Several future 
interferometers have plans to image the shadows around 
massive black holes: e.g. VSOP-2 for radio (Hirabayashi et 
al. 2005), MAXIM for soft X-ray (see MAXIM web page: 



http:/ /maxim. gsfc.nasa.gov/ ). While these planned interfer- 



ometers will image the massive black hole, direct imaging of 
the stellar-mass black holes are difficult because the appar- 
ent size of the black hole shadow is extremely small. For ex- 
ample, in the case of the black hole with 10 solar mass at the 
distance of 10 pc, the apparent size of the black hole shadow 
(~ 5 Schwarzschild radii) is about 0.1 micro- arcsecond. So, 
even with the future planned radio or X-ray interferometer, 
the direct mapping of the shadow of the stellar-mass black 
hole is practically impossible. However, if the stellar-mass 
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black hole in the binary system is occulted by its companion 
star, the eclipsing light curves give the direct information for 
the spatial distribution of the brightness patterns of the ac- 
cretion flows around the black hole. This is the primary mo- 
tivation of the present study. In terms of the studies about 
the strong-gravity region around the stellar-mass black hole, 
some observations in the field of X-ray spectroscopy pro- 
posed the values of the spin parameters of stellar-mass black 
holes (e.g. Zhang et al. 1997, Miller et al. 2002). In addition, 
recently the works in this field have achieved the significant 
advances for understanding the accretion flows in the strong- 
gravity regions around the black holes (e.g. Davis et al. 2006, 
Done & Gierlinski 2003, 2006, Gierlinski & Done 2004). 

In the present study, we study the method to investi- 
gate the physical information included in the strong-gravity 
region around the black hole by light-curve analysis for an 
accretion disc around a black hole occulted by its compan- 
ion star. The basic idea for the eclipsing light curve analysis 
was first proposed by Fukue (1987). Watarai, Takahashi & 
Fukue (2005) investigated the eclipsing light curves for the 
supercritical accretion flows which may be suggested as the 
accretion flow model for some X-ray black-hole binaries. In 
principle, light curves obtained at the time of eclipsing have 
information about the black hole if the region around the 
black hole is optically thin. However, the real observational 
data of the eclipsing light curves may contain the additional 
effects which are not sometimes simple. One of such effects 
is the atmospheric effects of the companion star. Recently, 
Pietsch et al. (2006) give the observational data of the eclips- 
ing light curves for the X-ray binary M33 X-7 detected by 
Chandra. In soft X-ray, the atmospheric effects in the com- 
panion star may vanish the relativistic effects in the eclipsing 
light curves. In other words, the eclipsing light curves also 
contain the information of the atmosphere of the compan- 
ion star. Thus, the eclipsing light curve analysis will be used 
for the investigation of the atmosphere of the companion 
star. In the present study, we also investigate the compan- 
ion star's atmospheric effects in the eclipsing light curves of 
the black-hole binaries, and compare these effects with the 
effects of the black hole's rotation. 

In §2, we give the calculation method of the eclipsing 
light curves. The adopted accretion disk model, the assumed 
spectral model and the calculation method of the eclipsing 
light curves without the atmospheric effects of the compan- 
ion star are given in §2.1. The assumed atmospheric model 
of the companion star and the calculation method of the 
eclipsing light curves with the atmospheric effects of the 
companion star are given in §2.2. In §3, we calculate eclips- 
ing light curves without the atmospheric effects of the com- 
panion star. After we calculate the eclipsing light curves 
around rotating black holes in §3.1, the analysis by using 
the normalized light curves and the statistical quantities of 
the skewness and the kurtosis are performed in §3.2 and 
§3.3, respectively. In §4, we calculate eclipsing light curves 
with the atmospheric effects of the companion star. After we 
show the sample images of the black hole shadows with the 
atmospheric effects in §4.1 and the normalized light curves 
in §4.2, the eclipsing light curves are analysed with respect 
to the observed photon energy in §4.3, the black hole spin 
in §4.4, and the inclination angle between the rotation axis 
of the disc and the direction of the observer in §4.5. Related 
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Figure 1. Bolometric flux distribution (erg cm -2 s _1 ) for black- 
hole spins of a/M = (left panels) and 1 (right panels), and 
inclination angles of i = 60° (top panels) and 80° (bottom pan- 
els). The relativistic standard discs are adopted. The inner edge 
of the accretion disc is assumed to be the marginally stable cir- 
cular orbit. For non-rotating black holes, indirect images can be 
seen in the region within the inner stable circular orbit. The flux 
distributions are asymmetric due to the disc rotation and the 
frame-dragging effects. The shapes of the shadows casted by the 
black holes are asymmetric due to frame-dragging effects around 
the rotating black holes. 



topics are discussed in §5. We give concluding remarks in 
the last section, §6. 



2 CALCULATION METHOD FOR ECLIPSING 
LIGHT CURVES 

2.1 Eclipsing Light Curves without Atmospheric 
Effects 

In order to calculate a light curve during an eclipse, some 
assumptions were required. In the present study, we assume 
the relativistic standard disc (Novikov & Thorne 1973, Page 
& Thorne 1974) at equatorial plane. Based on this model, 
the radial profile of the effective temperature of the accre- 
tion flow can be calculated. Then, at each radius, we assume 
the black body spectrum with the effective temperature cal- 
culated from the relativistic standard accretion disc model. 

The calculation method for the observed spectrum of 
the emission in the Kerr metric is established (e.g. Cun- 
ningham & Bardeen 1973). At each light ray, the observed 
specific intensity, /„ obB [erg s -1 cm -2 str -1 Hz -1 ], of the ac- 
cretion disc observed at the photon frequency v ^ s is calcu- 
lated as I„ ohs = g 3 Iv rost , where u TCBt is the photon frequency 
measured in the local rest frame of the accretion flow, T e ff is 
the temperature of the relativistic standard accretion disc at 
the radius where the light ray intersects, i„ roat is the specific 
intensity measured at the local rest frame of the accretion 
disc, and g is the frequency ratio defined g = ^obs/^rest- 
The frequency ratio g is calculated same as Cunningham 
& Bardeen (1973). In this study, we assume i„ rost is the 
black body spectrum with the effective temperature at the 
radius where the light ray intersects, i.e. I Vlost = -B„ rest (T) , 
where T is the temperature of the accretion disc. By us- 
ing the relation di/ b s = gdv ras t and the assumption of the 
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(a) Ingress 




(b) Egress 




Figure 2. Schematic diagrams for the eclipsing black- hole bina- 
ries when the ingress (top panel) and the egress (bottom panel). 

isotropic radiation at the local rest frame of the accretion 
disc, the observed flux F b s [erg s _1 cm -2 ] is calculated as 
-Fobs = (V rf2 ) / g!v obB dv ohs = (1/d 2 ) / g i I Vimt dv rast where 
d is the distance to the accretion disc. By integrating the 
observed flux in x-y plane, we can calculate the observed 
bolometric luminosity L D oi [erg s ]. On the other hand, 
the observed flux at the observed photon frequency v b s , 
i^obL [erg s _1 cm -2 Hz -1 ], is calculated by using the delta- 
function as F„ obs = (1/d 2 ) J g 4 Iv IQBt S(v ohs -gv Icst )dv Icst . By 
integrating the flux in x-y plane, the luminosity at the ob- 
served frequency u b B , L Vohs [erg s -1 Hz -1 ] is calculated. 
The luminosity are calculated by the ray-tracing method, 
which is commonly used in the past studies (e.g. Fukue & 
Yokoyama 1988, Takahashi 2004). When we transform the 
luminosity from the local rest frame to the Boyer-Lindquist 
frame by calculations of tetrads, the relativistic effects are 
included, that is, the special relativistic effects such as beam- 
ing effects and Doppler boosting, and the general relativistic 
effects such as gravitational redshift, gravitational lensing 
and frame-dragging effects. The frequency ratio g includes 
these special and general relativistic effects. 

For all calculations in the present study, we assume a 
black hole mass of WMq and mass accretion rate of lAfEdd, 
and the disc thickness is neglected. Here MEdd is the Ed- 
dington mass accretion rate defined as MEdd = -t-Edd/c 2 = 
1.4 x 10 17 (M/M Q ) [g/s] where L E dd is the Eddington lumi- 



Table 1. Cross sections [cm 2 ] for photoionization absorption due 
to HI, Hel, Hell and Compton scattering in X-ray of 0.1 keV, 1 
keV and 10 kcV. 
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nosity (e.g. Kato, Fukue & Mineshige 1998). Also, we im- 
plicitly assume that the binary system has a relatively high 
inclination angle, i > 60°. In figure [1] we show examples 
of the spatial distribution of the observed flux of the accre- 
tion flows in the vicinity of the black holes. For all panels 
of figure [1] the mass center of the black hole is at the cen- 
ter of the image. The width and the height of the images 
in figure [T] are both 20 M where M is the gravitational ra- 
dius and 1 M corresponds to 0.5 Schwarzschild radius. We 
define z-axis and y-axis as the horizontal line and the ver- 
tical line crossing at the center of the each image of figure 
[U In figure [U the position of the mass center of the black 
hole is at (x, y) = (0, 0), and the displayed region is —10 
[M]< x < 10 [M] and -10 [M]< y < 10 [M] . Based on the 
calculations of the observed luminosity described above, we 
calculate the eclipsing light curves. 

In figure [2l we show the schematic diagrams for the 
eclipsing black- hole binaries when the ingress (top panel) 
and the egress (bottom panel). We also plot z-axis and y-axis 
in panels of figure [2] When the ingress and the egress phases 
described in figure [2] the companion star is crossing the line 
of sight between the observer and the black hole binary along 
the direction of x-axis. Since the size of the companion star 
is generally much larger than the emission region around the 
black hole, we assume the carve of the stellar surface as the 
straight line in our calculations. Although in our calculations 
the brightness profiles of the accretion disc are calculated 
within r < 200 [M], the light curves are calculated in the 
range of —60 [M]< x < 60 [M]. That is, in terms of the range 
in y-axis direction, we use the brightness profile of the disc 
in the range of -200 [M]< y < 200 [M], but the brightness 
profiles of the accretion disc are given within r < 200 [M]. 

2.2 Atmospheric Model of Companion Star and 
Eclipsing Light Curves with Atmospheric 
Effects 

In the present study, we use the atmospheric model calcu- 
lated for the solar atmosphere (Daw et al. 1995, Fontenla et 
al. 1993, Gabriel 1976). In the top panel of figure [3] we plot 
HI number density JVhi [cm -3 ] (solid line) and temperature 
T [K] (dashed line) of the adopted atmosphere model. Here, 
H [km] denote the height above the photosphere. Temper- 
ature of the atmosphere suddenly increase and the number 
density decrease in the transition zone between the lower 
chromosphere (~ 10 3-4 [K]) and the upper corona (~ 10 6 
[K]). The two separate temperatures of HI correspond to 
two stable regions of the cooling function of HI (e.g. chapter 
6 in Spitzer 1978). 

In figure [4] we show the schematic diagrams for the ob- 
server, the eclipsing black-hole binaries, the companion star 
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with its atmosphere and the light ray penetrating the atmo- 
sphere. Here, the light ray penetrates the atmosphere with 
the length of 2^ ma x. Umax is the length from the center of the 
companion star to the maximum height of the atmosphere, 
and H is the length between the center of the companion 
star and the point where the light ray penetrate the atmo- 
sphere. The optical depth, t(u), in the atmosphere can be 
calculated as 

t{v) = J dl (cr abs iV abs + cr sca iV sca ), (1) 

where <r a bs and er sca are cross sections of the absorption and 
scattering, respectively, and iV a b s and iV aca are the number 
densities of the absorber and the scattering medium, respec- 
tively. The integrations are performed along the atmosphere 
of the companion star in the line of sight. 

In soft X-ray band, the photoionization absorption by 
hydrogen and helium contribute to opacity. The photoion- 
ization absorption cross sections of HI, Hel and Hell in 
soft X-ray are summarized in table [T] We also calculate the 
Compton scattering cross section in the same table. Since 
for 0.1 keV and 1 keV the scattering cross section is sig- 
nificantly lower than the absorption cross sections, we can 
neglect the effects of the Compton scatterings. On the other 
hand, for 10 keV the scattering effects are dominated. In 
this case, we only take into account the out-going photons 
as the scattering effects for simplicity. We see that the pho- 
toionization absorption cross sections of HI, Hel and Hell 
are comparable in soft X-ray. Calculations by Fontenla et 
al (1993) shows that the number density of Hell, Nnell, is 
significantly lower, e.g. by order 4, than that of HI, Nm, 
in the solar atmosphere, i.e. -/VhgII <C Nm, and the number 
density of Hel is lower by about order 1 than that of HI 
through the atmosphere, i.e. NhcI ~ Nm. Therefore, in the 
present study, we neglect the scattering effects in the atmo- 
sphere for 0.1 keV and the absorption by Hell for 0.1 keV, 
1 keV and 10 keV. Then, we calculate the optical depth in 
the atmosphere as 

t{v) = 2 [ di (amNm + OHdiVHci + crcompAy, (2) 
Jo 

where £ m&K = (#i ax - H 2 fl 2 , # max is set to 10 5 [km], 
the number density of electrons N e is calculated only in the 
photosphere based on numerical table given by Fontenla et 
al. (1993). However, according to Fontenla et al. (1993) the 
number density of electrons in the atmosphere is lower than 
that of HI by factor 3-4 in the almost region of photosphere, 
then the scattering effects are negligible in the photosphere. 
Near the transition zone, the number density of the electron 
is lower than that of HI within factor 1, and for cases of 10 
keV the scattering effects contribute to the optical depth. 
By using the optical depth calculated above, we solve the 
radiative transfer equation along the atmosphere in the line 
of sight. The brightness decreases along the light rays in the 
atmosphere by the exponential of the optical depth. 

By using the optical depth r calculated above, the ob- 
served luminosity with the atmospheric effects is calculated 
by multiplying the factor e -T ^ obB ' to the observed lumi- 
nosity without the effects of the atmospheric effects. The 
observed luminosity without the atmospheric effects corre- 
sponds to the case of r = 0. Since the distance between the 
accretion disc and the atmosphere of the companion star is 
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Figure 3. HI number density TVjii [cm -3 ] (solid line) and tem- 
perature T [K] (dashed line) of the adopted atmosphere model 
(top panel). The atmosphere model has been constructed for the 
solar atmosphere (Daw et al. 1995). The number density of Hel, 
Nn e i, is assumed to be Nn e i = O.IJVhi- The corresponding trans- 
parency ratio, e~ r , of the modeled atmosphere (bottom panel) for 
X-ray: 0.1 keV (long dashed line), 1 keV (solid line), 10 keV (short 
dashed line). Here, r is the optical depth of the photoionization 
absorption by HI and Hell. 

very large, here we evaluate the optical depth r at the ob- 
served photon frequency v h s - In the bottom panel of figure 
we show the corresponding transparency ratio, e~ T , of the 
modeled atmosphere for X-ray: 0.1 keV (long dashed line), 1 
keV (solid line), 10 keV (short dashed line). In the cases of 1 
keV and 10 keV, while the lines in the corona are optically 
thin, the photoionization absorption is effective in the pho- 
tosphere. By using the optical depths calculated as figure [3] 
we calculate the eclipsing light curves with the atmospheric 
effects of the companion star in later section. 

In the present study, we do not include contributions 
from the heavy elements which are important in some cases 
because the photoelectric cross sections are generally pro- 
portional to where Z is the proton number. However, we 
do not intend to give the quantitatively precise information 
of the eclipsing light curves in this section, but to roughly 
investigate the possibility whether the atmospheric absorp- 
tion effects can be negligible or not. Moreover, even if we 
calculate more precise light curves including the effects of 
the heavy elements, there is no observational targets having 
the solar-type atmosphere for which the precise light curves 
can be used to fitting the observational data (See also dis- 
cussion for the case of M33 X-7). 
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Figure 4. Schematic diagrams for the observer, the eclipsing 
black-hole binaries, the companion star with its atmosphere and 
the light ray penetrating the atmosphere. Here, the light ray pen- 
etrate the atmosphere with the distance 2£ max - J^max is the dis- 
tance from the center of the companion star to the maximum 
height of the atmosphere, and H is the distance between the cen- 
ter of the companion star and the point where the light ray pen- 
etrate the atmosphere. 
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Figure 5. Eclipsed light curves before (left panel) and after (right 
panel) a total eclipse by a companion stars. For each panel, the 
lower abscissa shows the time of eclipse and the upper abscissa 
shows the position of eclipse with respect to the position of black 
hole (x = 0). The inclination angles are 60° (solid lines) and 
80° (dashed lines). The black hole spins are 0, 0.5, 0.9 and 1 from 
bottom to top for i = 60° and 80° , respectively. The lines showing 
zero time (time= s) indicate the boundary of a complete eclipse. 
We assumed an orbital velocity of 200 km s _1 , a black hole mass 
of 1OM0 and mass accretion rate of IM^d^. 



3 ECLIPSING LIGHT CURVES WITHOUT 

ATMOSPHERIC EFFECTS OF COMPANION 
STAR 

3.1 Calculated Light Curves 

By using the observed bolometric luminosity Lboi, in figure 
[5]we show the light curves at the ingress (left panel) and 
egress (right panel) of an eclipse. The observed bolometric 
luminosity is the function of the position of the surface x s of 
the companion star, i.e. Lboi = iboi(is). The observed bolo- 
metric luminosity is calculated by integrating the observed 
flux in x-y plane. For the ingress phase with the stellar sur- 
face at x, the bolometric luminosity Lboi (a?) is given as 



L bol (x) = 



200[M] /■200[M] 

dx' I dy' 

200[M] 



du T , 



4 T 

9 Iv 



(3) 



where g and I„ rost are the functions of the light ray, x and 
y', and we set the observer's distance to be d — 100 [M]. As 
denoted in the previous section, we assume I Viest = B Vlost . 



In the same way, for egress phase with the stellar surface at 
x, 

r W0[M] 



Lboi(x) = — / dx' I dy' I du [cst g 4 h^t - (4) 

J-200[M] J-200[M] J 

We calculate the observed bolometric luminosity for the 
position of the stellar surface x in the range of — 60 [M] 
<x<60[M}. 

The observed bolometric luminosity can be also the 
function of time. Here, we define the origin of the time when 
the stellar surface is at 60 [M] for the ingress and at —60 
[M] for the egress. The region between ±60 [M] are enough 
to cover the effectively emitting region of the accretion disc 
observed at 1 keV, so we use these values in this study when 
we calculate the normalized light curves. In figure [S] the up- 
per and lower abscissas are the position of the stellar surface 
x and the time t. The larger values of the bolometric lumi- 
nosity is achieved for the larger spin parameter of a black 
hole. A rotating black hole transfers angular momentum to 
the accreting gas via the frame-dragging effect which is also 
calculated when the frame transformation by using tetrads. 
This effects increases the observed luminosity. 

Also, the marginally stable circular orbit r ms decreases 
as the black hole spin increases, e.g. r ms = 6 M for a/M — 
and r ms = M for a/M = 1 (Bardeen, Press & Teukolsky 
1972). Then, the black hole rotation increase the observed 
luminosity and the emission region in the vicinity of the 
black hole. The angular velocity uo of the zero angular mo- 
mentum observer (ZAMO) with respect to distant observer 
is calculated as u) = —gt<t>l 'g<j>4> = s'^/ff" = 2Mar/A. Here, 
A = (r 2 + a 2 ) 2 ~ a 2 A sin 2 9 where A = r 2 - 2Mr + a 2 . In 
this study, we assume 6 = 90° for the photon emitting re- 
gion. Since this angular velocity uj is roughly proportional 
tow oc 2M«/r 3 , the effects of the black hole rotation dras- 
tically decrease for large r. 



3.2 Normalized Light-Curve Analysis 

Although in figure we clearly see the effects of the black 
hole rotation as the increase of the bolometric luminosity, 
it is known that the mass accretion rate also increase the 
bolometric luminosity for some inclination angles (Watarai, 
Takahashi & Fukue 2005). So, in addition to the maximum 
bolometric luminosity, we analyze the shapes of the light 
curves which contain the information of a black hole spin. 
In order to do this, we define the normalized bolometric 
luminosity as L bo i(x) / Lho\(x = ±60M). That is, for the 
ingress phase, 



Lho\{x) 
L bo i(-60) 



200 (-200 [■ 

dx I dy' \ ripest 5 4 h xeBt 

J -200 J 
200 (-200 r 

dx' I dy' I ripest g 4 h taat ,(5) 

60 J -200 J 

where the unit for the length is M, and for the egress phase, 
Lboi (a;) 



Lboi(+60) 



dx' / dy' / digest g Iir TBat 

200 J -200 J 

60 /.200 <• 

dx' / dy' / df rC st g 4 L/ re8 t(6) 

-200 J-200 J 

Since we calculate the eclipsing light curves in the region 
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Figure 6. Normalized eclipsed light curves before (left panels) 
and after (right panels) a total eclipse by a companion stars. For 
each panel, the lower abscissa shows the time of eclipse and the 
upper abscissa shows the position of eclipse with respect to the 
position of black hole (x = 0). The ordinate shows the normal- 
ized bolometric luminosities defined as L-^ \(x) / ' L-^ Q \(x = 60M) 
for the ingress and L^ i(x) / L^, i(x = — 60M) for the egress. Here, 
L bo i(x) is the observed bolometric luminosity when the stellar 
surface is at x. That is, for the ingress phase, we normalize the 
observed bolometric luminosity L^ \(x) by the observed bolomet- 
ric luminosity when the stellar surface is at L\y \(x = 60M), and 
for the egress phase, the observed bolometric luminosity is nor- 
malized by the observed bolometric luminosity when the stellar 
surface is at L-t, \(x = 60M). In the same way, the observed lu- 
minosity L v (x) [erg s — 1 Hz -1 ] can be also normalized by the 
observed luminosity when the stellar surfaces are at x = ±60M. 
The inclination angles are are 60° (top panels) and 80° (bottom 
panels), and the black-hole spins are a/M = (solid lines), 0.5 
(long dashed lines), 0.9 (short dashed lines) and 1 (dotted lines). 
The lines showing zero time (time= s) indicate the boundary of 
a complete eclipse. We assumed an orbital velocity of 200 km s _1 , 
a black hole mass of 10A/q and mass accretion rate of lAfrjdd- 



of -60M < r < 60M in the present study, L b ol(a: = T60M) 
represent the maximum bolometric luminosity at the ingress 
and the egress, respectively. Figure [6] shows the calculated 
normalized light curves at the ingress (left panel) and egress 
(right panel) of an eclipse. From this figures, we see that the 
black hole rotation clearly deforms the shapes of the nor- 
malized light curves. For larger spin parameters, the light 
curves changes in the closer region of the black hole center 
(x — 0[A/]). This is because the effective emission region 
of photon from the accretion disc become more compact for 
larger spin parameter as shown in figure[T] Then, the drasti- 
cally changed timescale of the normalized luminosity become 
shorter for the larger spin parameter. We also investigate the 
effects of the inclination angle i. At first glance of the upper 
panels and the lower panels of figure [6] there is little depen- 
dence on the inclination angle. In order to see the effects of 



Figure 7. Same as figure [6] The inclination angles are are 60° 
(solid lines), 70° (long dashed lines), and 80° (short dashed lines) 
and 85° (dotted lines). The black-hole spins are a/M = (top 
panels) and 1 (bottom panels). 

the inclination angles on the shapes of the light curves, in 
figure [7] we plot the normalized light curves for the black 
hole spin of a/M = (upper panels) and 1 (lower panels). 
Although the absolute value of the bolometric luminosity 
depends on the inclination angle largely as we can see in fig- 
ure[5l from figure[7]the effects of the inclination angle on the 
normalized light curves is very small. Then, the shape and 
the variation timescale of the normalized light curve are the 
possible candidates of the good indicator of the black hole 
spin. 

3.3 Skewness and Kurtosis 

In addition to the normalized light curves given in the pre- 
vious section, we present the other indicator of the effects 
of the black hole spin. We introduce the statistical quan- 
tities of skewness, S, and kurtosis, K, as indicators of the 
light-curve asymmetry. These are the same as the analysis 
in Watarai, Takahashi & Fukue (2005). The skewness and 
kurtosis represent the deviation of the observational data 
from the Gaussian (normal) distribution. Skewness, S, is a 
measure of the degree of asymmetry of a distribution. On 
the other hand, kurtosis, K, is the degree of peaky feature 
of a distribution. 

We have done the skewness and kurtosis analysis for 
the inverted bolometric light curves as calculated in figure 
From the data of the inverted normalized light curve, we 
first obtain as /, = Lbol.max — iboi (ti) where U is the time 
of each mesh. The statistical quantities are defined as 

71 

p* = h/J2 fi > (7) 

!=() 
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where n is the number of data points, Pi is the normalized 
distribution calculated from the bolometric light curves as 
shown in figure [5] i is an weighted average time and a is the 
standard deviation. If the left tail is more pronounced than 
the right tail, the skewness becomes negative (S < 0). If 
the reverse is true, the skewness becomes positive (S > 0). 
The skewness for a normal distribution is zero S = 0, and 
any symmetric data should have a skewness near zero. The 
kurtosis for a standard normal distribution is three. Posi- 
tive kurtosis indicates a peaked distribution and negative 
kurtosis indicates a flat distribution. 

In figure [Sj we show the dependence of the calculated 
skewness, S (left panel) and kurtosis, K (right panel) on 
the black hole spin, a/M, for eight inclination angles (60°, 
65°, 70°, 75°, 80°, 82°, 85° and 88°). Since the marginally 
stable orbit become larger for non-rotating black hole than 
the rotating cases, then the emission region become wider 
than the rotating one. Then, for larger spin parameters, the 
skewness becomes smaller value. This is clearly seen in the 
left panel of figure [S] On the other hand, for larger spin 
parameters, the emission region and the variation timescale 
of the light curves become smaller. Then the light curves are 
expected to be more peaky than cases of the non-rotating 
black hole. Then, the kurtosis become larger value for the 
larger spin parameters of black holes. This feature can be 
seen in the right panel of figure [5] While the kurtosis cannot 
be good indicator of the mass accretion rate for non-rotating 
black holes as shown in Watarai, Takahashi & Fukue (2005), 
both of the skewness and the kurtosis can be the possible 
indicator of the black hole spin. 



4 ECLIPSING LIGHT CURVES WITH 

ATMOSPHERIC EFFECTS OF COMPANION 
STAR 

While in the last section we do not take into account the 
atmospheric effects of the companion star, in this section, 
we calculate the eclipsing light curves with the atmospheric 
effects. One of the key questions in terms of the atmospheric 
effects are whether the information of the black hole spin in 
the eclipsing light curves is completely smeared out by the 
atmospheric effects. One of the main goals in this section is 
to answer this question. We also investigate how the eclips- 
ing light curves depend on the observed photon energy in 
X-ray and the inclination angle. 



4.1 Black Hole Shadows with the Atmospheric 
Effects 

We first calculate the images of the black hole shadows with 
absorption effects in the atmosphere of the companion star. 
In figure [9] we show images of black hole shadows when 
ingress at the observed frequency of 1 keV for cases with no 
absorption effects (top panels) and atmospheric photoioniza- 
tion absorptions (bottom panels). We assume the mass of the 
companion star to be 1 Mq. The size of each image is same 
as figure Q] i.e. -20 [M]< x <20 [M] and -20 [M]< y <20 
[A/]. For no absorption case for top panels, the edges of the 
companion star are located at -20, -10, and 10 in units 
of M (left to right). For the cases with the atmospheric ef- 
fects (bottom panels), we define these edges as positions with 
t(v) — 7. We can see that the brightness around the black 
hole shadow are decreased by the photoionization absorp- 
tion in the atmosphere. In addition to the decrease of the 
absolute values of the brightness, the clear edges of the com- 
panion star which can be seen for no absorption cases are 
smeared out by the absorption effects. This is because the 
effective region with absorption effects is comparable to the 
size of X-ray emitting regions. 

4.2 Normalized Light Curves with the 
Atmospheric Effects 

We next calculate the eclipsing light curves with effects of 
the photoionization absorption. In figure 1101 we show nor- 
malized light curves with (thick lines) and without (thin 
lines) the atmospheric effects of the companion star. The 
observed energy is set to be 1 keV in figure 1101 The black 
hole spins are a/M = (solid lines) and 1 (dotted lines), and 
the observed inclination angles with respect to the rotation 
axis of the accretion disc are i — 60° (top panels) and 80° 
(bottom panels). While the normalization light curves with 
no absorption effects largely change in a short time when 
the most luminous parts of the accretion disc are occulted 
or unocculted by the sharp edge of the companion star, the 
normalization light curves with absorption effects gradually 
decrease or increase compared to the cases with no absorp- 
tion effects. This is because the atmosphere of the compan- 
ion star absorb some fraction of X-ray photons gradually 
before the most luminous part of the accretion disc is com- 
pletely occulted or unocculted. In other words, the width of 
the absorption region which causes the gradual changes of 
the X-ray brightness of the accretion disc is comparable to 
size of the luminous parts of the accretion disc as denoted 
above. 

The differences between the normalized light curves 
with absorption effects for a/M = and 1 become smaller 
than the cases of no absorption effects. When egress, i.e. 
t > 0, the normalized light curves with a/M — rise more 
rapidly than the cases of a/M = 1. That is to say, the 
timescale of the increase of normalized L„ for a/M — 
is longer than that for a/M. This is because the most lumi- 
nous parts of the accretion disc for a/M = 1 is more compact 
than those for a/M. This results in the shorter crossing time 
of the companion star along the luminous parts of the disc 
for a/M = 1 than the cases for a/M = 0. This feature can 
be seen in the case of no absorption by the companion star. 

Since the size of the most luminous parts of the ac- 
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Figure 8. Skewness and kurtosis for different black-hole spins and inclination angles, i = 60°, 65°, 70°, 75°, 80°, 82°, 85° and 
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Figure 9. Images of black hole shadows with a/M = 0.5 when ingress at the observed frequency of 1 keV for cases with no absorption 
effects (top panels) and atmospheric absorptions (bottom panels). The size of each image is same as figure [TJ i.e. —20 < x/M < 20 and 
—20 < y/M < 20. For no absorption case for top panels, the edges of the companion star are located at -20, -10, and 10 in units of M 
(left to right). For the absorption case for bottom panels, we define these edges as positions with t(u) = 7. 



cretion disc depends on the black hole spin, the crossing 
timescale of the companion star along the most luminous 
parts of the accretion disc can be one of the key physical 
quantities which distinguish the cases of rapidly rotating 
black holes and the cases of no rotating black hole. But these 
effects are generally smeared out by the atmospheric effects 
of the companion star at some fraction. In our calculations 
assuming the solar atmosphere, the atmospheric absorption 
can not completely smeared out the difference between the 
cases of rapidly rotating black holes and the cases of no ro- 
tating black hole. Then, the observed eclipsing light curves 
at 1 keV with the atmospheric effects of the companion star 



possibly contain the information of the black hole spin. How- 
ever, since the effects of the atmosphere is larger than the 
effects of the spin, only when we have the reliable atmo- 
spheric model, we can determine the black hole spin from 
the eclipsing light curves observed by the wavelength where 
the atmospheric absorption is effective. 



4.3 Observed Photon Energy 

Here, we investigate the observed energy dependence of the 
atmospheric absorption effects in the images of the black 
hole shadows and the normalized light curves. In figure [TT1 
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Figure 10. Normalized light curves with (thick lines) and with- 
out (thin lines) the atmospheric effects of the companion star at 
observed energy of 1 keV for black hole spins a/M = (solid 
lines) and 1 (dotted lines), and the observed inclination angles 
with respect to the rotation axis of the accretion disc i = 60° 
(top panels) and 80° (bottom panels). 



Figure 12. The normalized light curves with effects of the atmo- 
spheric absorption for observed energy of 0.1 keV (dashed lines), 
1 keV (solid lines) and 10 keV (dotted lines). The black hole spin 
is a/M = and the inclination angles of the observers are i = 60° 
(top panel) and 80° (bottom panel). 



images of the black hole shadows in the accretion discs are 
shown for observed photon energy of 0.1 keV (left panel), 1 
keV (middle panel) and 10 keV (right panel). Here, we do not 
include the atmospheric absorption effects. In these panels, 
the inclination angle of the observer is i — 80° and the black 
hole spin is a/M = 1. The absolute values of the luminosi- 
ties are largest for the case of 1 keV. On the other hand, the 
larger the observed photon energy becomes, the more com- 
pact the regions of the luminous parts in the accretion discs 
become. Especially, in the case of 10 keV, the size of the 
luminous parts in the accretion disc is smallest and the lu- 
minous parts exist in the blue shifted parts of the accretion 
disc which are located in just outside the black hole shadow. 
As Takahashi (2004) calculated, the size and the positions 
of the black hole shadows in the accretion discs clearly de- 
pend on the black hole spins. Thus, the eclipsing light curves 
for higher photon energy can contain the information of the 
black hole shadows in case that the atmospheric smearing 
effects are negligible. 

The normalized light curves with effects of the atmo- 
spheric absorption are shown in figure[l2]when a/M = and 
in figure [T3l when a/M = 1 for observed energy of 0.1 keV 
(dashed lines), 1 keV (solid lines) and 10 keV (dotted lines). 
The black hole spin is a/M — and the inclination angles 
of the observers are i = 60° (top panels) and 80° (bottom 
panels). The differences between cases of 0.1 keV, 1 keV and 
10 keV are clearly shown in these figures. The higher the ob- 
served photon energy become, the shorter the timescales of 
the changes of the normalization light curves become. This 
is because the larger the observed photon energy becomes, 
the more compact the regions of the luminous parts in the 



accretion discs become and then the crossing time of the 
companion star is shortest in the case of 10 keV and longest 
in the case of 0.1 keV. These features can be seen in both 
cases of a/M — in figure [12] and 1 in figure [13] 



4.4 Black Hole Spin 

In order to see the effects of the black hole spin in the nor- 
malized light curves when the observed photon energies are 
0.1 keV and 10 keV, we plot the normalized light curves 
same as figure [TU] but for the observed photon energy of 0.1 
keV in figure ITU and 10 keV in figure IT51 In the case of 0.1 
keV in figure [141 the difference of the normalized light curves 
between the cases of a/M = and the cases of a/M — 1 
are very small for both cases of i = 60° and 80°. Thus, 
it seems very difficult to see the information of the black 
hole spin in the observed light curves. On the other hand, 
in the case of 10 keV, we can clearly see the difference of 
the normalized light curves between the cases of a/M = 
and the cases of a/M = 1 for both cases of i — 60° and 80°. 
The timescale of the change of the normalized light curve 
for the case of a/M — 1 is longer than the timescale for the 
case of a/M = for both cases of ingress and egress. The 
normalized light curves are gradually decrease and increase 
because of the atmospheric effects and the most luminous 
parts of the accretion discs exist in just outside of the black 
hole shadows as denoted above. Moreover, the size of the 
black hole shadows become smaller for larger spin and the 
position of the black hole shadow are shifted to the opposite 
side of the most luminous parts of the accretion discs (Taka- 
hashi 2004). The timescale when the most luminous parts 
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Figure 11. Images of the black hole shadows in the accretion discs for observed energy of 0.1 keV (left panel), 1 keV (middle panel) 
and 10 keV (right panel). In these panels, the inclination angle of the observer is i = 80° and the black hole spin is a/M = 1. 
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Figure 13. The same as figure H"2l but for a/M = 1. 

can be seen by the observer become longer when a/M = 1 
than the case of a/M = because the most luminous parts 
exist in more inner regions when a/M = 1 than the case of 
a/M = 0. 



4.5 Inclination Angle between the Rotation Axis 
and the Direction of the Observer 

We also investigate the dependence of the normalized light 
curves on the observed inclination angles, i, with respect to 
the rotation axis of the accretion disc. While so far we plot 
the normalized light curves for both cases of i = 60° and 
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Figure 14. The same as figure 1101 but for the observed photon 
energy of 0.1 keV. 



80°, in order to see the dependence on the inclination an- 
gle i, we plot the normalized light curves with same black 
hole spins in the same panel. The normalized light curves 
for cases of a/M = (top panel) and 1 (bottom panel) are 
shown in figure [T6l For each panel, the lines are plotted for 
the combinations of the inclination angle i = 60°, 70° and 
80° and the observed photon energy to be 0.1 keV, 1 keV 
and 10 keV. While the observed images of the black hole 
shadows in the accretion discs are clearly different for the 
different inclination angles i, the dependences on the incli- 
nation angles are very small in the normalized light curves. 
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Figure 15. The same as figure llOl but for the observed photon 
energy of 10 keV. 
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5 DISCUSSION 

5.1 Eclipsing Light Curves and Observational 
Feasibility 

Recently, Pietsch et al. (2006) give the clear observational 
data of the eclipsing light curves for the X-ray binary M33 
X-7 detected by Chandra. While the absorption effects in the 
companion atmosphere seem to vanish the relativistic effects 
in the eclipsing light curves in their observational data, our 
calculations show that the information of the black hole spin 
in the eclipsing light curves can not be completely smeared 
out for the observed photon energy of 1 keV and 10 keV if we 
assume the solar-type atmospheric model. However, as de- 
noted above, even if we have the reliable observational data 
with high signal-to-noise ratio, only when we have the reli- 
able atmospheric model of the companion star, the physical 
information of the strong-gravity regime, such as, a black 
hole spin, can be obtained. 

In addition, in the present study we use the solar atmo- 
spheric model as a basic sample of the atmospheric struc- 
ture. In the case of the eclipsing X-ray binary M33 X-7, the 
companion star is an 06III star. Although some recent ob- 
servational data suggest the existence of corona in the atmo- 
sphere of Herbig Ae/Be stars (Stelzer et al. 2006), it remains 
one of the unsolved issues whether the high temperature 
stars such as O-type stars have corona in their atmosphere. 
If there is the corona in the atmosphere of O-type stars, 
the stellar atmosphere will have two-temperature structure 
such as the solar atmosphere. In such cases, our calculations 
roughly calculate the light curves for the case of M33 X-7. 

Basically, the count rates of the X-ray photons from 
X-ray binaries in external galaxies with the possible excep- 
tion of the Magellanic Clouds are generally much too low to 
provide the precise eclipse mapping of the inner portions of 
the accretion discs even in Chandra observations. In addi- 
tion, the time resolutions for the data of the eclipsing light 
curves for M33 X-7 are insufficient to precisely map the small 
strong- gravity regime in the center of the accretion flows. 
Even in the case that the calculations in the present study 
are not used to the present low count rate observations in X- 
ray, future X-ray data with higher count rates obtained by 
the future X-ray telescope such as Constellation-X Obser- 
vatory in NASA (see, http://constellation.gsfc.nasa.gov/) 
or XEUS (X-Ray Evolving Universe Spectrometer) in ESA 
(see, http://www.esa.int/science/xeus) will provide the in- 
formation of the extremely strong-gravity regions of the ac- 
cretion discs. 



5.2 Accretion Disc Model 



Figure 16. The normalized light curves for cases of a/M = 
(top panel) and 1 (bottom panel). For each panel, the lines are 
plotted for the combinations of the inclination angle i = 60° , 70° 
and 80° and the observed photon energy to be 0.1 keV, 1 keV 
and 10 keV. 



While our calculation assumes the standard accretion disc, 
for the accretion discs with larger mass accretion rate, 
Watarai, Takahashi & Fukue (2005) pointed out the light 
curves are deformed by the self-occultation effects, and 
shows that for larger mass accretion rate the skewness be- 
comes small value and the kurtosis is nearly constant. On 
the other hand, our analysis shows the larger the black hole 
spin becomes, the smaller the skewness becomes and the 
larger the kurtosis becomes for any observed inclination an- 
gles. The dependence of the kurtosis on the black hole spins 
contrast with the dependence of the mass accretion rate. 
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We expect that the effects due to the black hole rotation 
can not be seen for discs with large mass accretion rate be- 
cause of the self-occultation effect. The same analysis for the 
accretion discs with high mass accretion rate around rotat- 
ing black hole is a topic for future study. Moreover, from 
the eclipsing light curves, it may be possible to examine the 
disk structure such as stress at the marginally stable orbit, 
the angular dependence of the emitted radiation, and the 
dynamical structure relating the magnetic processes such 
as magnetorotational instabilities in the supersonic region 
which are the fundamental unsolved issues about the ac- 
cretion flow models. These issues are also topics for future 
studies. 



5.3 X-ray Heating Effects 

So far, we implicitly neglect the effects of the heating of the 
atmosphere due to the absorbed X-ray photons in the atmo- 
sphere of the companion star. According to the calculations 
of the structure of the heated atmosphere (e.g. Basko & 
Sunyaev 1973), the temperature of the atmosphere increase 
by about factor 2 due to the absorbed X-ray photons. The 
number density is also changed by similar factor. Although 
these changes causes the quantitative changes of the optical 
depth in X-ray of the atmosphere of the companion star, we 
do not expect the qualitative changes of the dependence of 
the eclipsing light curves on the physical parameters. Since 
the absorption effects basically do not change the timescale 
of the change of the eclipsing light curves, the effects of the 
X-ray heating will not drastically alter the dependence of 
the crossing timescale of the luminous parts of the accretion 
disc. In addition to the changes of the optical depths, a X- 
ray heating effects causes a radiatively driven stellar wind 
originated from the X-ray heated atmosphere (e.g. Blondin 
1994, Wojdowski et al. 2001). The atmospheric profile of the 
companion star is not the same between ingress and egress 
due to this effect. Since this effects also change the optical 
depth of the atmosphere in X-ray, we expect the effect will 
not change the timescale of the change of the eclipsing light 
curves drastically. 

5.4 Cases for Neutron Stars and Seyfert Galaxies 

Asymmetric brightness distribution exists not only in black 
hole X-ray binaries, but also in neutron star X-ray binaries 
(NSXBs). This is because the relativistic effects still works in 
NSXBs. Even in the case that the emission from the neutron 
star is dominant, the asymmetric component still remains in 
its light curves as long as the emission from the neutron star 
surface is isotropic. The observational data for the eclipsing 
neutron star binaries are actually obtained (e.g. Homan et 
al. 2003). In the case of optically thin along the light paths, 
we can estimate the emitting radius from the duration time 
of ingress or egress, and measure the spin of the neutron star 
from the skewness analysis of the asymmetric eclipsing light 
curves. According to our results, the result does not depend 
on the inclination angle. Weaver & Yaqoob (1998) shows 
the "deep minimum" in the X-ray light curves of Seyfert 1 
galaxy MGC-6-30-15. Some obscuring body in the line of 
sight seem to cause the occultation phenomena in the light 
curves. Only when such events occur many times and the 



typical features of the light curves are obtained, and the 
photons from the vicinity of the black hole are directly de- 
tected, it may be possible to obtain the physical information 
in the vicinity of the black hole, such as the size of emission 
region and black hole spin. 



6 CONCLUDING REMARKS 

In this paper, we propose an eclipsing light-curve diagno- 
sis for accretion flows around rotating black holes. When 
emission from an inner part of the accretion disc around 
black hole is occulted by a companion star, the light curves 
at ingress and egress show the asymmetric features due to 
the relativistic effects such as Doppler boosting originally 
pointed out by Fukue (1987). Our calculations include the 
asymmetric feature due to the black hole rotation. After 
we give the calculation method of the eclipsing light curves 
with and without the atmospheric effects of the companion 
star in §2.1 and §2.2, respectively, the numerical calculation 
of the eclipsing light curves are performed and the analysis 
for the curves are given. When no atmospheric effects of the 
companion star, the effects of the black hole rotations on the 
light curves are clearly seen in the light curves (§3.1), the 
normalized light curves (§3.2) and the statistical quantities 
of skewness and kurtosis (§3.3). We show that the eclipsing 
light curves without the atmospheric effects clearly reflect 
the effects of the black hole's rotation. We also investigate 
the eclipsing light curves with the atmospheric effects of the 
companion star. Based on the adopted atmospheric model, 
after evaluating the absorption and scattering effects in X- 
ray, we calculate the images of the black hole shadows in 
§4.1, the normalized light curves §4.2. The analysis of the 
eclipsing light curves with the atmospheric effects are given 
in terms of the observed photon energy (§4.3), the black hole 
spin (§4.4) and the inclination angle (§4.5). The higher the 
observed photon energy become, the shorter the timescales 
of the changes of the normalization light curves become. This 
is because the crossing time of the companion star is short- 
est in the case of 10 keV and longest in the case of 0.1 keV. 
Since the size of the most luminous parts of the accretion 
disc depends on the black hole spin, the crossing timescale of 
the companion star along the most luminous parts of the ac- 
cretion disc can be one of the key physical quantities which 
distinguish the cases of rapidly rotating black holes and the 
cases of no rotating black hole. In our calculations assuming 
the solar atmosphere, while the effects of black hole spins 
are negligible when the observed photon energy of 0.1 keV 
for both cases with no atmospheric effects and with atmo- 
spheric effects, for the observed photon energies of 1 keV 
and 10 keV the atmospheric absorption can not completely 
smeared out the difference between the cases of rapidly ro- 
tating black holes and the cases of no rotating black hole. 
However, these effects are smeared out by the atmospheric 
effects of the companion star at some fraction, and in our 
atmospheric model, the effects of the atmosphere are much 
larger than the effects of the black-hole spin. Therefore, even 
in the case that the light-curves contain the information of 
the black hole spin, it may be difficult to extract the infor- 
mation of the black hole spin if we do not have the realistic 
atmospheric profiles, such as, the temperature, the number 
densities for several elements. Even in such cases, based on 
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our calculations, the light curve asymmetries due to the ro- 
tation of the accretion disc exist. Only when we have the 
reliable atmospheric model, in principle, the black hole spin 
can be determined from the eclipsing light curves observed 
around 1 keV. 
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